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Abstract 

The escape rate of a stochastic dynamical system can be found as an expansion in powers of the noise strength. In 
previous work the coefficients of such an expansion for a one-dimensional map were fitted to a general form containing 
a few parameters. These parameters were found to be related to the fractal structure of the repeller of the system. 
The parameter a, the "noise dimension", remains to be interpreted. This report presents new data for a showing 
that the relation to the dimensions is more complicated than predicted in earlier work and oscillates as a function of 
the map parameter, in contrast to other dimension-like quantities. 
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1. Introduction 

Physical systems all experience noise to some de- 
gree. Whether from unknown degrees of freedom, 
thermal fluctuations, quantum fluctuations or many 
other causes the dynamics are affected in all systems 
from neuronal systems jl|2j to cosmology [3] . Differ- 
ential equations governing physical systems are ide- 
alised and in reality will be affected at some scale. 
This motivates the study of systems which evolve 
with noise, specifically maps with additive noise. 

The problem of whether a particle escapes or at 
what rate particles escape from an open system 
is another subject that has many physical appli- 
cations, from an asteroid escaping a planet 'A\ in 
celestial mechanics to the chemical reactions that 
occur in transition state theory [5]. The study of 
billiards has expanded to include open systems |6l7j 
and open optical billiards may have industrial ap- 
plications [8 . All of these physical systems may also 
experience noise. 

This paper considers asymptotic expansions for 



the escape rate from one dimensional unimodal 
maps with small additive noise. One-dimcnsional 
maps have shown to be models of properties such 
as universality |9ll0j . intermittency |11I12| and de- 
terministic diffusion |13I14| . The effect of noise on 
the properties of one-dimensional maps has been 
well studied |15ll6ll7ll8ll9l20j particularly in the 
vicinity of bifurcation and crisis points [21122123] . 
Properties such as escape rates and lifetimes of 
transient behaviours have been shown to have a 
nonmonotonic dependence on noise strength |24| . 
In contrast to most of the previous work we discuss 
the high order terms in the weak noise expansion 
which determine the rate of convergence or diver- 
gence and hence the usefulness of the expansion. 

The escape rate from a one-dimensional map with 
additive noise can be calculated as an expansion in 
terms of the noise strength using information about 
the periodic orbits of the corresponding determinis- 
tic map. Different approaches exist [25126] but the 
most accurate results are obtained using a numerical 
method described in [57]. This method accurately 
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yields a large number of terms in the expansion but 
takes a large computational time. In [28j the ex- 
pansion is reformulated into a new asymptotic form 
where the parameters are related to properties of 
the system including the complex fractal dimension. 
Many physical systems show complex dimensions. 
Ref. [in] discusses the structure of the bronchial tree 
in the mammalian lung and Ref. 30] reviews a num- 
ber of other applications including the rate of escape 
from a stable attractor surrounded by an unstable 
limit cycle as a function of noise strength and the 
formation of black holes with mass smaller than the 
Chandrasekhar limit. 

However, the parameter a, the "noise dimension" , 
remains to be found. Finding the significance of this 
parameter could enable the expansion of the escape 
rate to be calculated more easily or to higher order. 
Conversely, it could also allow the properties of a 
system to be inferred from experimental data for the 
escape rate. 

In |28j it is suggested that the noise dimension 
may be related to the generalised Renyi dimensions 
of the fractal repeller which are known to change as 
the map changes. This paper uses numerical tech- 
niques to calculate the noise dimension for a range 
of maps and compares the data to the Renyi dimen- 
sions in an attempt to find the significance of the 
remaining parameter. Generalisations to higher di- 
mension or continuous time systems are likely, how- 
ever, this will require a greater understanding of one- 
dimensional systems. 

Section 2 reviews some results from periodic orbit 
theory which lead us to the escape rate for determin- 
istic and stochastically perturbed systems. In Sec- 
tion 3 the theory of asymptotic expansions is gener- 
alised to include complex "exponents" and the new 
formulation of the expansion of the escape rate is in- 
troduced. Numerical results for the noise dimension 
are found and discussed in Section 4 and finally our 
conclusions arc given in Section 5. 

2. Periodic orbit expansions for the escape 
rate 

The systems considered include unstable and 
chaotic trajectories, particularly when there is noise, 
so we study the evolution of probability densities 
rather than individual trajectories. The evolution 
of a density of particles moving in a determinis- 
tic dynamical system is governed by the Perron- 
Frobenius evolution operator |31j . If the dynamics 



are described by a discrete time map Xn+i = f{xn) 
where Xfi IS Oil real number, then the probability 
density, 4>{x), evolves as 

0«+i(2/) = (-C o 0„)(y) = J (j)n{x)C{y,x)dx (1) 

with the evolution operator C given by 

C{y,x) = Siy-f{x)). (2) 

The Perron-Frobenius evolution operator can also 
be generalised to higher-dimensional and continuous 
systems. 

For a sufhcicntly uniformly chaotic system with a 
uniform distribution of initial conditions the number 
remaining in the considered region at long time n 
is proportional to e"'''", where 7 is the escape rate. 
The escape rate of an open deterministic system can 
easily be calculated from the leading eigenvalue, ly, 
of the evolution operator using the relation, 7 = 
— In ly. Other properties such as averages, Lyapunov 
exponents and correlation functions are also related 
to the leading eigenvalue of the evolution operator. 
This eigenvalue is simply the leading root of the 
spectral determinant of the linear operator which 
can be expressed as an expansion in terms of powers 
of traces of the evolution operator as shown below. 

0= det(l - zC) = exptrln(l - zC) 
= l-ztrC-—itrC^~{trCf ) + ■■■. (3) 

Here z is the inverse eigenvalue and 

trC" ^ J dxC'{x,x) = J dx5{x~ r{x)). (4) 

The only values of x which contribute to this integral 
are those that lie on periodic orbits. Periodic orbit 
theory provides a simple formula to calculate these 
traces using information about the periodic points 
of the system |32|33j . 

In the case of a stochastic dynamical system where 
the system experiences some additive noise the 
equivalent evolution operator is the Fokker-Planck 
evolution operator |32j . If the dynamics are de- 
scribed by a discrete time map Xn+i = f{xn) + cr^n 
where a is the noise strength and ^„ are instances 
of a random variable with normalized distribution 
p(^) then the probability density, (f){x), evolves as 
Eq. ([1]) but now the evolution operator £ is given by 

C{y,x) = 5„{y-f{x)) (5) 



2 



where 



Jim 



(6) 



As in the deterministic case, generahsations to 
higher-dimensional and continuous time systems 
are possible. Equivalent results are shown in [34j 
where the trace formula is found for a noisy flow. 

For stochastic systems the problem of finding the 
escape rate can also be reduced to finding the traces. 
Here we have 



5a{xi - f{xQ))...S^{xQ - f{Xn-l)) 

xdxo---dxn~i 



(7) 



but now there is no simple formula that can be 
used from periodic orbit theory. Two analytical 
approaches have used Feynman diagrams [23] and 
smooth conjugations [IB] to express the traces as 
expansions in terms of the noise strength. However, 
they have only obtained results up to order cr^. In- 
stead a different more numerical approach may be 
taken which expresses the evolution operator as a 
matrix which is calculated in the vicinity of each 
periodic orbit to find the traces and then obtains 
the eigenvalue as an expansion in terms of the noise 
strength [57] . When used with high precision arith- 
metic this method can obtain results of up to order 
64 in the noise strength [28] . With such a high order 
expansion obtained to a high degree of accuracy it 
is possible to look at the asymptotic form of the 
late terms. Reformulating the expansion in terms 
of new parameters could allow it to be calculated 
more efficiently or to a higher order. 

Previous work in [28] proposed one such new for- 
mulation. By studying one particular map most of 
the new parameters were found but a parameter 
known as a remained. However a tentative relation 
between a and the dimension was proposed. The 
aim of this paper is to present new data for a by 
looking at different maps and to show that the be- 
haviour is more complex than predicted in the pre- 
vious work. A brief description of the derivation of 
the parameters follows but a more thorough expla- 
nation is given in the original paper [55] . 

3. The expansion for the escape rate in 
powers of the noise strength 

The expansion for a smooth map with symmetric 
noise is fitted in [28] to the following asymptotic 
form, 



m 

where 

^ oo 

X! Cn{m + a + in(3)l 



V2 



(8) 



(9) 



Here a is the noise strength and V2m are coeffi- 
cients that can be obtained using the local matrix 
approach [27]. 

Eq. ^ without the complex inf3 (and hence sum 
over n) is familiar from the theory of asymptotic 
series for integrals of the form / exp((7(u) /a)h{u)du 
[35136] . The quantity 



_2m 



(10) 



would normally be expected to tend to a constant 
as m tends to infinity, however, in [28] it was found 
to exhibit oscillatory behaviour. The complex in/3 
were added in order to replicate this behaviour. 

The parameter ctq, called the singulant, is given 
by the variation of g{u) between the point defin- 
ing the asymptotic expansion and the closest criti- 
cal point [35^. Here this corresponds to finding the 
probability of returning to the repeller from the crit- 
ical point of the map [28] . The integration variable 
u is in our case the whole path {xi\ and a maxi- 
mally probable path is then a critical point in this 
infinite-dimensional space; see Section 4 for the ex- 
plicit calculation. 

The remaining parameters, a, [3 and c„, were 
found in [55] to be related to the fractal structure of 
the set of repelling points of the map. We motivate 
this connection as follows: in [37| a similar analytic 
expansion was found for the electrostatic potential 
near the edge of a middle-third Cantor set with a 
uniform charge on it. The middle-third Cantor set 
contains two copies of itself scaled down by a factor 
of three so has dimension d = In 2/ In 3. The value 
of the potential a distance ^ from the edge of the set 
is given by a series of where the d values are 
the complex solutions of the equation that describes 
the dimension of the set, 3*^ = 2. 



d = 



In 2 27rm 



In 3 In 3 



(11) 



where n is an integer, a is the dimension and (3 is re- 
lated to the spatial scaling factor of the fractal. Os- 
cillatory terms arise in this expansion as a result of 
the complex values of d corresponding to n 7^ 0. In 
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this expansion the coefhcients, c„, have a rapid ex- 
ponential decay so only \n\ < 2 are required. Com- 
plex fractal dimensions such as these appear in many 
physical applications |29|30j . 

Looking back at the original expansion defined by 
dl]) and ([U we can analogously relate the values of a 
and (3 to the properties of the fractal repeller of the 
map. The middle-third Cantor set is an exactly self- 
similar fractal, however, the repeller of the map may 
not be exactly self-similar and may be a multifractal. 
In this case the spatial scaling factor would vary 
with position. In [5H] f3 is shown to be related to 
the spatial scaling factor at the most probable point 
of return to the considered region from the critical 
point and is given by 

where / is the map and is the point of return. 
However, for the parameter a, the noise dimension, 
it is not so simple. As the repeller may be a multi- 
fractal, we need the concept of generalised dimen- 
sion. The Renyi dimension of a fractal covered by 
boxes of size e is defined as 



^hmi^ 

q — 1 e^o Ine 



(13) 



where pi is the amount of the fractal in box i [32] . 
The middle-third Cantor set is exactly self-similar 
so has dimension d ~ In2/ln3 for all values of q. 
However, when the repeller for the map is not ex- 
actly self-similar there are infinitely many different 
dimensions. In this context it is not clear which di- 
mension is related to a. 

4. Results: ot vs. Renyi dimensions 

Recall from the previous section that the escape 
rate of a stochastically perturbed map may be fitted 
to the form given by ([8]) and ([9]) . Explicit expressions 
can be found for the parameters ctq a-nd (3. However, 
the parameters a and c„ remain to be found. In this 
section we calculate a and c„ numerically by fitting 
the new form to the initial expansion and minimis- 
ing the least squares error after finding the values 
of the other parameters. Following the equation for 
the potential near the middle-third Cantor set we 
assume that the c„ values decay exponentially and 
also use only \n\ < 2 which is found to be justi- 
fied when the first c„ values are calculated. We then 
compare the data to the dimensions in an attempt 
to better understand the significance of a. 



We have considered systems where the noise is 
additive and the noise strength is given by cr, so an 
individual trajectory will evolve with time as 



(14) 



where the one-dimensional, real variable x„ is the 
position ait — n, f is the map and ^„ are instances of 
an independent and identically distributed random 
variable. Previous work [27l28j has concentrated on 
the map where 



/(^) = -(l-(l 
a 



2x) 



(15) 



and a = 0.8. The map can be changed by vary- 
ing the map defining parameter, a, between and 
1. Fig. [T] shows this map for a = 0.8. These spe- 




Fig. 1. The map defined by l|15|l when a = 0.8. 

cific maps were chosen because it is simple to find 
the inverse branches and they have complete binary 
symbolic dynamics which means periodic orbits are 
easy to enumerate using inverse iteration, however, 
the method does not require complete symbolic dy- 
namics. With this map we are interested in the [0, 1] 
region and in order to escape a particle must leave 
this region. The noise chosen here is Gaussian with 
distribution 



V2^ '' 



(16) 



however, any other peaked smooth symmetric noise 
distribution would be appropriate including those 
that are space dependent. 

For small a the probability of returning to the 
repeller from the critical point of the map where 
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f'{x) = is given by (ignoring algebraic depen- 
dence) 



exp 



(17) 



where we identify ctq with the singulant that appears 
in ([9]) ^28j. In this hmit the probabihty is given to 
leading order by the probability of the most hkely 
trajectory that returns to the repeller from the crit- 
ical point. For the map defined by the critical 
point is Xc — 0.5. From here the most likely point 
of return to the repeller is at the edge of the set. If 
the trajectory returns after one step = 1 or if it 
returns after more than one step x,. = 0. The prob- 
ability of a trajectory {xq, xi, x„} is given by 



exp 



20-2 



(18) 



Finding the most probable trajectory with xq = 
Xc = 0.5 and equivalent to minimising 

the sum over all possible lengths n and values of 
xi,X2-, ...jXri-i- Calculating this minimum for each 
n shows that longer trajectories are more likely. Pre- 
vious work in |28j found that when a = 0.8 in (|15p 
the probability is maximised by the infinite trajec- 
tory {0.5, 1.00244613635157, -0.00024587488150, 
2.460023246 x 10-^ -2.46015082 x 10^^^ 
- 2.4601636 X 10-^-2.460165 x 10-^...,0} and 
CTo = 32.31850341240166-1/2 Qn this trajectory 
the approach to the repeller is geometric with ratio 
1/10 = l/f'{xr) which is the inverse of the spatial 
scaling factor. 

Using the most likely point of return, a;,. = 0, we 
can also find (3 as follows. 



2tt 



27r 



27r 



\nf'{xr) ln(f(l-2x,)3) ln(8/a)- 



(19) 



In [55] one specific case of the map was considered 
where a — 0.8. When the fitting described above 
was carried out a value of a = —1.290 was found. 
The Renyi generalised dimensions [32] were also cal- 
culated for this system and Z?oo was found to be 
0.2957. At this point a tentative connection between 
a and Doc was proposed, a = —1 — Doo- This re- 
port presents further data for a, as such conclusions 
cannot be drawn from only one piece of data. 

Results obtained are shown in Figs. [2][5l Fig. [2] 
shows a plot of the behaviour of a as the map is 
varied along with —1 — Doo- The data for a was ob- 
tained using a noise expansion with terms up to or- 
der 48 calculated using periodic orbits with length 



up to 12. The new form is fitted to the initial expan- 
sion using coefficients of cr^^ and higher order terms 
to perform a least squares fit. 




Fig. 2. Plot of the noise dimension, a, and —1 — Doo for the 
map defined by I I15I I. 
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a 

Fig. 3. Plot of Doo, -Do and D-oo for the map defined by 

Comparison between a and —1 — Doo shows that 
the graphs are clearly different to each other and the 
relation proposed in [28j holds only approximately 
for 0.8 < a < 0.9. 

As a tends to zero the repeller of the map tends 
to that of an infinitely high tent map. The repeller 
of the tent map is a Cantor set, which is an exactly 
self-similar fractal. All the dimensions will have the 
same value and this value will tend to zero as the 
height increases. As a result of this we would expect 
the plot of a to tend to a constant as a tends to 
zero if it were a simple function of dimension, but 
instead, we find it oscillates increasingly fast as it 
approaches zero. 

The variation of the generalised Renyi dimensions 
Doo, Dq and I?_oo is shown in Fig. (3] The plot of any 
other Renyi dimension should also be smooth and 
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lie between the plots of Doo and I?-oo- Examination 
of this plot shows that the dimensions do tend to 
the behaviour of the tent map as a tends to zero, 
however, they do not exhibit behaviour as oscillatory 
as that of a and so it seems unlikely that there is a 
simple relationship between a and Dq for any q. 

In the case of the exactly self-similar fractal all 
of the Renyi dimensions were identical and a had 
a simple relationship to the dimension but in the 
case studied here the fractal repeller is not exactly 
self-similar so the Renyi dimensions are all different. 
In this situation it is not clear which dimensions 
to use and this may account for the more complex 
behaviour. 

We can also apply this method to different maps. 
Fig.[4]shows data for a selection of different maps of 
the form 



f{x)^-{l-{l-2xf). 
a 



(20) 



As the order of the map increases a increases and 
the oscillations become larger. Fig. [5] shows the cor- 
responding variation of the dimension Dq for these 
maps. It is clear that the variation of the dimensions 
with a or fc is very different to the variation of a. 
This further supports the idea that a is not a simple 
function of the dimensions. The only clear similarity 
is that as k increases the plots get closer together, 
however, this is expected as the map changes less as 
k increases. 




Do ° 




Fig. 4. Plot of the noise dimension, a, for maps defined by 
l(20)l. 
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a 

Fig. 5. Plot of Do for maps defined by I I20I I. 

fitted this expansion to a new form with parame- 
ters that depend on the fractal properties of the re- 
peller. In this work we have focused on one param- 
eter, the noise dimension a, which was thought to 
be related to the complex fractal dimension [29l30j 
but remained to be found. 

Although we have not found the noise dimension, 
a, explicitly in terms of other known quantities we 
have produced new data showing how it varies as 
the map varies. We have calculated a as the map 
parameter is changed for the quartic map given by 
(fT5|) and higher order maps given by ([20)1 . We have 
shown that any relation to the Renyi dimensions 
of the repeller is much more complicated than pre- 
dicted. The dimensions are smooth monotonic func- 
tions of the map parameter a whereas a oscillates. 

Future work will continue to look for the signifi- 
cance of the noise dimension and its relation to the 
static and dynamic properties of the repeller. Find- 
ing this parameter could allow the expansion for the 
escape rate to be found more simply or allow us to 
draw conclusions about the properties of a system 
from experimental measurements of the escape rate. 
The data presented for a poses many new questions, 
such as, what is the limiting behaviour as the height 
or power of the map tends to infinity and can the os- 
cillations in a be described as a function of the map 
defining parameter a? We hope that answering these 
questions will provide new insight into the nature of 
the escape from stochastic dynamical systems. 



5. Conclusions 

The escape rate from a stochastically perturbed 
one-dimensional map can be given as an expansion 
in terms of the noise strength. Previous work [28j 
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